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We present an extensive density functional theory study of the electronic, magnetic and transport 
properties of GaAs and AlAs digital ferromagnetic heterostructures. These can be obtained by 
5-doping with Mn the GaAs layers of a GaAs/ AlAs superlattice. Our analysis spans a range of Mn 
concentrations and considers the presence of compensating defects such as As antisites. 

In the defect-free case all the heterostructures studied present an half-metallic electronic structure. 
In contrast when As antisites are present the half-metallic state is destroyed and the heterostructures 
behave as dirty planar metals. In this case they show a large p-type metallic conductance in the 
Mn plane mainly due to majority spin electrons, and an n-type hopping-like conductance in the 
GaAs planes mainly due to minority spin electrons. This suggests that if the As antisites can 
be kept far from the Mn planes, spatial separation of the different spin currents can be achieved. 
Finally we show that in the case of AlAs/(Ga,Mn)As digital ferromagnetic heterostructures the 
AlAs/GaAs valence band offset produces an additional confining potential for the holes responsible 
for the ferromagnetism. Therefore the ferromagnetic coupling between the Mn ions becomes larger 
and more robust to the presence of As antisites. 

PACS numbers: 75.50.Pp, 71.20.Nr, 71. 15. Mb 



I. INTRODUCTION 

(Ga,Mn)As pj is the prototype of a new class of mag- 
netic materials named diluted magnetic semiconductors 
(DMS) 013. These are obtained by doping ordinary 
semiconductors with transition metals. In the case of 
(Ga,Mn)As the Mn ions occupy the Ga sites and provide 
both localized spins (5=5/2) and holes. The spin holes 
are then antiferromagnetically coupled to the Mn ions 
and this gives rise to hole-mediated long range ferromag- 
netism via a Zener-like mechanism 

The potential impact of this material on the semicon- 
ductor industry is huge, since it opens the possibility of 
adding ferromagnetism to the AlAs/GaAs system, an im- 
portant step toward the implementation of the spin de- 
gree of freedom in an electronic device |5J . To date several 
concept devices have been demonstrated, including spin 
polarized light emitters @ and electrically controlled fer- 
romagnetism |2J, and one can envision (Ga,Mn)As among 
the building blocks for a spin-based quantum computer 
scheme @- 

Despite these indisputable successes the critical tem- 
perature (T c ) of (Ga,Mn)As hardly exceeds 110 K 0,0 
and this poses severe limitations to future commercial 
applications. At present low temperature annealing pro- 
cessing has allowed an increase of the Curie temperature 
from 110 K to about 170 K 0, ^j, but room tempera- 
ture ferromagnetic (Ga,Mn)As appears difficult to pro- 
duce. However these recent annealing experiments have 
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demonstrated an important point: the low critical tem- 
perature of (Ga,Mn)As is not an intrinsic limitation of 
the material, but is affected by the inability to control 
the Mn and holes concentration independently. 

A possible way to improve the control over the elec- 
tronic and magnetic interactions is to produce magnetic 
semiconductor superlattices. These are the digital ferro- 
magnetic heterostructures (DFH) , which are obtained by 
(5-doping with Mn a low temperature GaAs MBE-grown 
layer [TlJ . Here Mn concentrations as high as 50% are ob- 
tained in a few (typically 2-3) GaAs monolayers and one 
may expect correspondingly higher large Curie temper- 
atures. However these structures show properties rather 
different from those of their random alloy counterparts. 
Here we report the most relevant experimental findings. 

1) The Curie temperature is rather low (~ 50 K) and 
independent of the separation between the Mn planes 

. T c usually decays with increasing GaAs layer thick- 
ness separating the MnAs sub-monolayers, and saturates 
for thicknesses larger than ~50 GaAs monolayers. The 
saturation is unexpected according to the mean field 
model for three dimensional systems, since the total Mn 
concentration in the sample decreases with the increase 
of the GaAs thickness jj| . This separation dependence 
suggests that DFH behave like planar systems. 

2) Hall measurements in the direction parallel to the 
MnAs planes show an anomalous Hall effect for undoped 
samples, which disappears upon Be doping [l2l [l3| . 
Large Shubnikov de Haas oscillations are found in doped 
samples, although surprisingly the carrier densities ex- 
tracted from the Hall coefficient and from the Shubnikov 
de Haas oscillations are different. This suggests that two 
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different carrier types could be present in the system. 

3) There is a correlation between metallicity and Curie 
temperature. In low Tq samples the transport is given 
by activated hole conduction and this is consistent with 
variable range hopping in two dimensions |l4i ITU Il6j| . In 
contrast, in the only DFH to date showing Tc around 
room temperature the transport is p-type and metal- 
lic ^tJ- However in these latter structures made from 
(Ga,Mn)Sb two phases may be present 0,0], with a 
diluted phase responsible for the ferromagnetism below 
40 K and a zincblende MnSb phase responsible for the 
room temperature ferromagnetism. 

4) GaAs/AlAs band engineering and spatial selec- 
tive doping [l8L fli^] allow the enhancement of the 
T c in AlAs/(Ga,Mn)As DFH with respect to their 
GaAs/(Ga,Mn)As counterparts. This enhancement is 
correlated with the enhanced hole concentration in the 
Mn layers and to the achievement of metallic conduc- 
tance. 

Finally we point out that some of the aspects described 
above are common to other 5-doped structures. For in- 
stance it has been recently demonstrated [2(j that Be 
5-doped low temperature GaAs undergoes an insulator 
to metal transition as the Be concentration is enhanced. 
In this case the transport changes from n-type thermally 
activated, to p-type metallic. The first is reported for 
small Be concentrations and is due to the hopping be- 
tween As antisite levels in the GaAs region, while the 
second dominates at large Be concentrations and is due 
to extended hole states in the Be rich region. 

From this brief overview it is clear that DFH present 
rather rich and complex physics, which calls for an exten- 
sive theoretical analysis. So far the magnetic properties 
of DFH have been studied only within the mean field ap- 
proach pll I22]] . while the transport has been investigated 
solely in the ballistic limit for the case of 100% Mn doping 
in the plane |23j| . Here we report an extensive ab-initio 
study of the electronic, magnetic and transport prop- 
erties of both GaAs/(Ga,Mn)As and AlAs/(Ga,Mn)As 
DFH, for various Mn concentrations and As antisite dop- 
ing levels. 

The paper is organized as follows. In the next section 
we describe our computational technique and we moti- 
vate the approximations made. Then we investigate the 
properties of GaAs/(Ga,Mn)As DFH, the effects of As 
antisites, and the properties of AlAs/(Ga,Mn)As DFH. 
Finally we conclude and we suggest new ways to manip- 
ulate the properties of DFH. 



II. COMPUTATIONAL TECHNIQUE 

We perform density functional theory (DFT) [24| cal- 
culations within the local spin-density approximation 
(LSDA). The use of LSD A for DMS is very well docu- 
mented [25|, and it provides a good description of the 
main physics of (Ga,Mn)As. Recently we have demon- 
strated [26} that self-interaction corrections to the LSDA 
do not strongly affect the band structure of (Ga,Mn)As, 
although they lead to strong localization and orbital or- 
dering of the Mn d shell in (Ga,Mn)N. For this reason we 
choose to work within the LSDA. 

Our numerical implementation, contained in the code 
SIESTA [27], 13 , uses pseudopotentials and a highly op- 
timized localized atomic orbital basis set. These two 
aspects make SIESTA extremely suitable for handling 
systems with a large number of atoms in the unit cell 
without a significant loss of accuracy. The drawback is 
that both the pseudopotentials and the basis set must be 
accurately optimized. 

First we consider the pseudopotentials. We use well- 
tested scalar relativistic Troullier-Martins pseudopoten- 
tials [23 with non linear core corrections [3(| and 
Kleinman-Bylander factorization [3]]]. The eigenvalues 
for the valence electrons of the free atom are compared 
with those generated for an all-electron calculation for 
different atomic and ionic configurations. Then we per- 
form total energy calculations for elementary solids com- 
paring the lattice constant, the bulk modulus and the 
band structure with reference calculations. These are 
performed with a well-converged basis set. Note that this 
is quite a delicate procedure, since with localized orbital 
basis sets the variational principle is not governed by a 
single-parameter such as the cut-off energy with plane- 
waves. 

Finally an optimized basis set is selected. The ba- 
sis functions in SIESTA are the product of an angular 
function with a given angular momentum, and a ra- 
dial numerical function. This latter is constructed as 
the DFT solution of the free pseudo-atom with an ad- 
ditional hard-wall confining potential. Furthermore, in 
order to enhance the variational freedom, several radial 
functions ("zetas") for the same angular momentum are 
constructed with the 'split valence' scheme (j^. In the 
case of (Ga,Mn)As the crucial aspect is to introduce sev- 
eral zetas for the Mn d orbitals, and the criterion we have 
adopted is that of reproducing the physics of the d shell 
of Mn in MnAs. More details are given in reference |33j . 
Here we only mention that the same procedure has been 
adopted for the Al pseudopotential and basis set. These 
have been checked for both metallic Al and zincblende 
AlAs. The reference configuration for the pseudopoten- 
tial is 3s 2 3p 1 with cut off radii 1.90, and 1.80 a.u. respec- 
tively for the s and p shells. Finally the basis set for Al 
has two basis functions for both s and p electrons, with 
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the same cut-off radius of 6.0 a.u. and the 'split norm' 
parameter is 0.15 [28j. 

All the calculations presented here are performed 
within a supercell scheme. Our supercell is constructed 
from a 2 x 2 x 3 zincblende cubic cell (lattice constant 
ao=5.65A) and contains 96 atoms in total. We mimic a 
DFH by replacing Ga with Mn ions only in one of the 
GaAs planes. We use periodic boundary conditions in all 
directions sampling 18 fc-points in the supercell Brillouin 
zone. This corresponds to a (Ga,Mn)As/GaAs super- 
lattice in which the (Ga,Mn)As planes are separated by 
6 GaAs monolayers (16.95A). In this supercell the Mn 
ions can occupy only eight possible positions in plane. 
These are arranged into two simple cubic lattices trans- 
lated with respect to each other along the diagonal of 
the xy plane. Since the exact positions of the Mn ions is 
rather important in determining the electronic structure, 
these are schematically presented in figure ^ 



lattices are occupied (namely Al, A2, A3, and Bl). Here 
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FIG. 1: The eight possible positions of the Mn ions in the 
supercell. 



III. (Ga,Mn)As DFH 

In this section we analyze the effects of the Mn concen- 
tration on the electronic properties of (Ga,Mn)As DFH, 
by calculating the band structure, the DFT total poten- 
tial, the strength of the FM coupling and the transport 
properties. Our main aim is to monitor the evolution 
of these quantities as a function of the Mn concentra- 
tion in the (Ga,Mn)As plane. In particular we want to 
establish whether there is a correlation between the Mn 
concentration and the metallicity of the system. 



A. Band structure 

In figure [21 we present the band structure of our 
(Ga,Mn)As/GaAs DFH for Mn concentrations of 12.5%, 
25% and 50%. In figures Et> and|2k, only the sites be- 
longing to one of the two cubic sub-lattices in the plane 
are occupied. This maximizes the mean Mn-Mn separa- 
tion. In contrast in figure |2Ji sites belonging to both the 
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FIG. 2: Band structure of the (Ga,Mn)As/GaAs DFH de- 
scribed in the text as a function of the Mn concentration. 
The directions are parallel to the MnAs planes. The upper 
panels correspond to the majority spin-band and the lower to 
the minority. The Mn concentrations are a) 12.5%, b) 25%, 
c) and d) 50%. In the case of 50% concentration we consider 
two different arrangement of the Mn atoms: c) Al, A2, A3, 
A4, d) Al, A2, A3, Bl. 

we plot the bands only along two directions parallel to 
the (Ga,Mn)As plane, since in the perpendicular direc- 
tion these are very similar to the case of MnAs planes 
embedded in GaAs |23j and they do not change signifi- 
cantly with the Mn concentration. In contrast for direc- 
tions parallel to the Mn planes, we expect a transition 
from the band structure of a (Ga,Mn)As random alloy to 
the metallic band structure of a zincblende MnAs plane 
|23[ . This is indeed the behavior observed in figure El If 
one increases the Mn concentration from 12.5% to 50% 
(from (a) to (b) to (c)) the Fermi level shifts downward 
in energy ending up deep in the majority spin valence 
band. At the same time also the spin-splitting of the va- 
lence band increases. Note that the DFH behaves as an 
half-metal at every concentration. The magnetic moment 
of the supercell therefore is always 4^b x Amh where Ny[ n 
is the number of Mn atoms in the cell. Mulliken popu- 
lation analysis [U 01 shows an orbital population for 
the Mn d shell of ~ 4.7 for the majority spin electrons 
and of ~ 0.8 for the minority, although both depend on 
the specific spatial arrangement of the Mn ions. These 
aspects are consistent with the picture of Mn as a single 
acceptor in GaAs. The Mn ions are in a d 5 state with an 
associated antiferromagnetically coupled hole, as for the 
case of the random alloys |3^| . 

However the situation of figure |2i is in stark contrast 
with this picture. In this case the Mn concentration is 
still 50% but three of the four Mn ions in the plane occupy 
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nearest neighbor positions. The band gap in the majority 
spin band closes and the material is an half metal with 
a completely metallic majority spin band. This suggests 
that the actual position of the Mn atoms in the plane 
is crucial in determining the electronic properties. The 
same sensitivity of the electronic structure to the posi- 
tion of the Mn ions is also present in the random alloys 
|35| . However this generally does not lead to a strong 
distortion of the band structure, and the bands obtained 
for Mn ions diluted in the supercell or occupying nearest 
neighbor positions are quite similar. In DFH the planar 
arrangement of the Mn ions makes the system more con- 
fined, and therefore more sensitive to inhomogeneities. 
This of course drastically affects the scattering properties 
of electrons in the (Ga,Mn)As planes, as we will show in 
the following sections. 
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B. DFT potential 

One fundamental question for understanding the 
physics of DFH is: "are the spin carriers confined in 
the (Ga,Mn)As plane or do they spread over the GaAs 
spacer?" . In order to answer this question it is useful to 
investigate the behavior of the total DFT potential along 
the superlattice direction z. This of course has the same 
periodicity as the atomic lattice. However we are not in- 
terested in the potential at the atomic scale, and instead 
we perform a "macroscopic average" 36] . The macro- 
scopic average is obtained by first taking a planar aver- 
age and then by averaging the result over the period of 
a GaAs monolayer along the superlattice direction. The 
resulting z-dependent potential is that felt by an electron 
with long wave-length at a density small enough not to 
perturb significantly the potential. Therefore this can be 
interpreted as the effective z-dependent potential in the 
spirit of the envelope function approximation. 

In figure |3| we plot the macroscopic average for the 
total DFT potential, the Hartree potential and the charge 
density distribution as a function of the position along 
the superlattice direction for different Mn concentrations. 
Notice first that there is indeed a confining potential in 
the Mn plane. The charge density is unevenly distributed 
along the superlattice and accumulates in the Mn plane. 
This is consistent with t he p icture of the Mn ions that we 
have given in the past 33], in which the majority spin 
hole is nearly bound to the Mn ion while the minority 
feels a much weaker potential. 

The z-dependent total potential is strongly spin- 
dependent. Generally it has a double well structure, with 
two potential minima, located in the (Ga,Mn)As plane 
and in the GaAs spacer respectively. These are sepa- 
rated by a potential barrier which grows if the Mn con- 
centration in the plane increases. We further investigate 
the nature of the confining potential by calculating its 



FIG. 3: Macroscopic average of a) the total DFT potential for 
majority spins, b) the total DFT potential for minority spins, 
c) the Hartree potential and the d) charge density distribution 
as a function of the position along the superlattice direction, 
for different Mn concentrations (x is the number of Mn ions 
in the plane). 



evolution upon increasing the Mn doping. For this pur- 
pose we define AJ and A| as the energy minima in in 
the (Ga,Mn)As plane and the GaAs spacer respectively 
measured with respect to the top of the energy barrier 
(see figure 0J. a is the spin index and f (|) indicates 
the majority (minority) spin electrons. We denote as A^ 
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FIG. 4: Definition of A? and Af . 

and A?? the same quantities for the Hartree potential. 

In figure [S] we present these quantities as a function of 
the number of Mn ions in the Mn plane (the maximum 
number is 8). First we notice that for minority spins both 
A| and A^ mimic closely the behavior of the Hartree po- 
tential (figure Et, and Et) ■ This means that the DFT 
potential for the minority spins is largely electrostatic 
with small contributions from the exchange part. In con- 
trast Aj for the majority electrons deviates strongly from 
the behavior of its electrostatic component, and the de- 
viation increases as the Mn concentration is enhanced. 
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FIG. 5: Af , A£, Af and A 2 as a function of the Mn concen- 
tration in the (Ga,Mn)As plane. 



If we now compare Ai with A2 for the two spin species 
we find a rather different behavior. For the majority 
spin electrons (figure [SJi) Ai > A2 at all Mn concentra- 
tions. Furthermore the difference Ai — A2 increases with 
increasing the Mn concentration. This means that the 
minimum in the (Ga,Mn)As plane is always the lower of 
the minima for majority electrons and it becomes more 
stable as the Mn concentration increases. 

In contrast, for the minority electrons there is not a 
striking difference between Ai and A2 meaning that they 
are less confined in the (Ga,Mn)As plane (see figurc[3l). 
It is also interesting to point out that A2 > Ai at low Mn 
concentration but there is a cross over for doping of the 
order of 50%. This is suggestive of the fact that the con- 
finement of minority electrons switches from the GaAs 
to the (Ga,Mn)As region upon the increase of Mn con- 
centration. Of course there are no empty minority spin 
states at the Fermi level and therefore in these structures 
the minority spins contribute little to the electron trans- 
port. However in the case of additional doping (for in- 
stance by intrinsic As antisites, Asc a ) an impurity band 
may form at the Fermi level, opening a transport chan- 
nel in the minority sub-band. The potential described 
here suggests that those electrons will hardly distribute 
over the (Ga,Mn)As plane at low doping but they will 
invade the (Ga,Mn)As region at large Mn doping. For 
this reason we expect a spatial separation of the spin 
currents (majority spin in the (Ga,Mn)As plane, and mi- 
nority spin in the GaAs spacer) for low concentrations, 
and strong spin-mixing at higher concentrations. 



C. Stability of the ferromagnetic state 

In this section we investigate the strength of the ferro- 
magnetic coupling in (Ga,Mn)As DFH, and in particular 
we discuss how the coupling depends on Mn doping. The 
relevant quantity to investigate is the energy difference 
Afa = Eaf — Efm between the total energy of the anti- 
ferromagnetic (£af) and the ferromagnetic (Efm) con- 
figurations of the supercell. Here we investigate the case 
of 25%, 37.5% and 50% Mn doping in the plane, where 
respectively two, three and four Mn ions are present in 
the supercell. This gives us some freedom to arrange the 
Mn ions in different ways and to investigate different an- 
tiferromagnetic configurations. A summary of the cases 
studied is presented in table [I] 

TABLE I: Summary table of the different magnetic configu- 
rations studied. The labels for the positions of the Mn ions 
are those of figure 

[ Mn] Configuration Position 



Af A2 A3 A4 Bf 



25% 


FMi 


T 








T 


25% 


AFMi 


T 








J 


25% 


FM 2 


T 


T 








25% 


AFM 2 


T 


I 








25% 


FM 3 


T 






T 




25% 


AFM 3 


T 






1 




37.5% 


FMi 


T 




T 


T 




37.5% 


AFMi 


T 




1 


1 




37.5% 


AFM 2 


T 




I 


f 




37.5% 


FMi„ 


T 




t 




T 


37.5% 


AFMl„ 


T 




T 




1 


50% 


FMi 


T 


T 


T 


T 




50% 


AFMi 


T 


I 


I 


T 




50% 


AFM 2 


T 


I 


T 


T 




50% 


AFM 3 


T 


T 


I 


1 




50% 


FM ln 


T 


T 


T 




T 


50% 


AFMl„ 


T 


T 


T 




I 


50% 


AFM 2n 


4 


T 


T 




T 


50% 


AFM 3 „ 


T 


T 


1 




1 



We ran a total energy calculation for each of the con- 
figurations of table n and fit our calculations to a third 
nearest neighbors Heisenberg model, in which the total 
energy E is expressed as 

nn 

E = - ^JijSi ■ Sj . (1) 

i<j 

Jij are the exchange constants, Si is the spin of the i-th 
Mn ions, and the sum runs up to third nearest neighbors. 
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The results of this fit are presented in table [HI where by 
convention we use \S\ = 5/2. Note that in our fitting 
procedure we have more configurations than parameters, 
and it is remarkable to observe that the fit is generally 
very good. From the table it is clear that the first nearest 

TABLE II: J\, J2 and J3 from the total energy calculations 
for the configurations of table [I] The units are meV. 



[Mn 


in plane 


Ji 


J2 J 3 




X 


= 0.25 


23.3 


2.9 5.6 


7.95 


X ■ 


= 0.375 


19.8 


1.4 4.9 


9.8 


X 


= 0.50 


13.3 


0.9 4.5 


9.35 



neighbor coupling constant, Ji, is responsible for most of 
the coupling which decays rapidly with the Mn-Mn sep- 
aration. It is also interesting to note that the second 
nearest neighbor coupling, J2, is small for all the concen- 
trations studied. Remarkably all the coupling constants 
are strongly dependent on the Mn concentration in the 
plane, and they decay quickly as this is increased. This 
decay, which is particularly severe for the first nearest 
neighbor coupling constant, is in stark contrast with the 
expectations of RKKY-like theories in which the J's in- 
crease as x 1 / 3 at zero temperature [32]] ■ However it is 
consistent with previous DFT-LDA calculations j3g for 
DMS random alloys, and illustrates once again the criti- 
cal dependence of the Mn-Mn ferromagnetic coupling on 
the local chemical environment. 

In order to put our calculations in perspective, next 
we estimate the ferromagnetic Curie temperature, Tq , of 
the system for different Mn concentrations. This can be 
simply obtained by using the mean field expression for a 
three dimensional Heisenberg model, which reads 

T c k B = ^S 2 n Mn ^2 Ji , (2) 

i 

where nun is the Mn concentration, and the sum runs 
over all the cation sites. This is a rather crude approx- 
imation and a complete thermodynamic theory should 
include the elementary spin-excitations |3^, |Zj|. More- 
over in the present case the Mn concentration is not a 
well-defined quantity since the (Ga,Mn)As region cannot 
be separated from the GaAs region. Roughly speaking 
one should consider the volume of (Ga,Mn)As to be the 
region around the Mn ions as thick as the range of the 
relative confining potential (see figure . This quan- 
tity is not clearly defined. Therefore, assuming that the 
(Ga,Mn)As volume does not depend on the Mn concen- 
tration, we prefer to evaluate only the following "mag- 
netic energy" 

^mag ~ % ^ ^ Ji ) (3) 
i 



which is proportional to the Curie temperature. 

From table[n]one notes that i? m ag has a non-monotonic 
dependence on the Mn concentration, presenting a max- 
imum for x = 0.375. Such behavior is generally observed 
in DMS random alloys, for which there is a maximum of 
Tc upon Mn doping, followed by a sharp decay for large 
Mn concentrations (above x = 0.05) This usually 

coincides with the loss of the metallic state. In contrast, 
in DFH made to date the situation seems to be reversed 
|TT| . with a larger Tc for larger Mn concentrations. This 
apparently contradicts our predictions. However two im- 
portant aspects need to be considered. First in actual 
DFH the transport is through variable range hopping 
|16| , while in our supercell calculations the system is "by 
definition" metallic. Secondly DFH usually present very 
strong compensation. This indicates that a large num- 
ber of donors, whose density is probably related to the 
in plane Mn concentration, are present in the DFH. In 
the next sections we will investigate systematically both 
the transport properties and the effects of the presence 
of donors. 



D. Ballistic Transport 

In this section we investigate the ballistic transport in 
DFH, with the aim of understanding the nature of the 
electronic states responsible for the conductance. The 
technique used is identical to that described in refer- 
ences |23t |41| and here we summarize only the main as- 
pects. The transport is calculated by using the self con- 
sistent tight-binding like Hamiltonian and overlap matrix 
computed by SIESTA. The matrix elements are obtained 
from the self-consistent charge density by evaluating nu- 
merically both two- and three-center integrals [2^] . Then 
we rewrite both the Hamiltonian and the overlap matrix 
in a tridiagonal form along the direction of the transport 
and we use periodic boundary conditions along the other 
directions. Finally, the fc-dependent transmission ma- 
trix t a {k) for the spin direction a is calculated by using 
our Green's function technique [4lT | , and the spin conduc- 
tance T a in the Landauer-Biittiker formalism [4^ . 

2 BZ 

r = ^Trt ff (fc) ( „Wt , (4) 

k 

where we integrate over the two-dimensional Brillouin 
zone in the plane orthogonal to the direction of the trans- 
port. Here we consider a two spin fluid model, where 
there is no mixing between the majority and minority 
spin currents. 

We study transport only in the direction parallel to the 
Mn plane. In the orthogonal direction in fact the trans- 
port is mainly due to hopping between the Mn planes 
and it is strongly suppressed if these are sufficiently far 
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apart [23. In figure we present the conductance as a 
function of energy for a 50% Mn supercell respectively in 
the FM state, with the Mn ions uniformly distributed in 
the plane (configuration FMi of table . 




E (eV) E (eV) 

FIG. 6: Conductance as a function of energy for a 50% Mn 
DFH in the FM state. The Mn ions are uniformly spaced in 
the plane (configuration FMi of table U} . The vertical line 
indicates the position of the Fermi level, (a) majority, (b) 
minority spins. 

As expected from the band structure (figure Et) the 
conductance as a function of energy shows an half- 
metallic behavior, with zero conductance for the minority 
spin band at the Fermi level. This 100% spin polariza- 
tion of the conductance persists down to 0.5 eV from the 
Fermi energy, where minority states at the top of the va- 
lence band start contributing to the current. Turning our 
attention to the orbital contribution to the conductance 
it is important to observe that, at the Fermi level the 
current is entirely due to a mixture of As p and Mn d ti 
states. It is also interesting to note that non-negligible 
Mn d contributions are present in the majority spin band 
for energies down to 4 eV below the Fermi level. This is 
an indication of the strong p-d hybridization in the ma- 
jority spin band, and in particular at its top. In contrast 
the Mn d contribution to the conductance is almost neg- 
ligible in the valence minority band, and Mn d states 
appear only for E > — 3 eV, relative to the bottom of the 
conduction band. 

This 100% spin-polarization at the Fermi level is very 
encouraging for the potential use of DFH as spin-injector 
for spintronics devices. However it is crucial to investi- 
gate how this feature survives when compensating defects 
are present in the system. 



IV. (Ga,Mn)As DFH: EFFECTS OF As Ga 

As in the case with the (Ga,Mn)As random alloys, 
DFH are also usually strongly compensated so that the 



hole concentration is considerably lower than the Mn con- 
centration. It is generally accepted that the strong com- 
pensation is due to donors, most likely of intrinsic defects. 
In particular both As antisites (AsGa) and inter- 
stitial Mn (Mn,) |43Ll44| have been indicated as relevant 
compensating defects. The relative abundance of those 
donors probably depends on the growth conditions, the 
Mn concentration, and the post-growth processing. Since 
DFH are usually grown under large As overpressure 
we believe that in this case As antisites dominate. Here 
we investigate how the electronic properties of a 50% Mn 
DFH changes upon Asca doping. 



A. Band Structure 

As in the previous section the band structure provides 
important information on the electronic properties of the 
DFH. In figure [7| we present the band structure for 50% 
in plane Mn DFH, where a single AsQ a is introduced 
into the GaAs spacer at midway between two consecu- 
tive (Ga,Mn)As planes (the total Asc a concentration is 
-2%). 




FIG. 7: Band structure for 50% Mn DFH where a single 
AsGa is introduced in the GaAs spacer at midway between 
two (Ga,Mn)As planes. The directions are: 1) Y-F-S in the 
(Ga,Mn)As plane, 2) T-Z perpendicular to the (Ga,Mn)As 
plane. The horizontal line indicates the position of the Fermi 
level. On the left (right) we plot the majority (minority) 
band. 

The most important feature of this band structure is 
that the half-metallic state is destroyed. This is due to 
the appearance of the Asc a level in the minority spin 
band, and to the fact that the Fermi energy is pined at 
this level. An As antisite in GaAs is a double donor with 
a doubly occupied deep level at midgap and an empty 
resonant state at the edge of the conduction band. These 



states have respectively the A and T2 symmetries of the 
Td point group 0. in (Ga,Mn)As the density of such 
defects is generally rather large and they can give rise 
to the formation of narrow impurity bands. In addition 
in the case of DFH it is likely that the As antisites will 
concentrate in the proximity of the Mn layers. This is 
confirmed by our total energy calculations, which show 
that there is an energy gain of approximately 0.9 eV when 
an As antisite moves from the middle of the cell to the 
Mn plane. 

The total magnetic moment of such a unit cell is about 
17.2/ie and depends weakly on the position of the As an- 
tisite with respect to the Mn plane (it is 17.33/ib when 
the Asca lies in the (Ga,Mn)As plane). Assuming a rigid 
band model and considering that, in absence of As anti- 
sites the magnetic moment of the cell is 16/ie, we con- 
clude that an Asc a contributes 0.4 electrons to the minor- 
ity band and 1.6 to the majority (we recall that one Asca 
also introduces an impurity state that can accommodate 
two electrons) . Therefore the presence of one As antisite 
in such a unit cell has two main effects: i) it compensates 
up to 1.6 holes in the majority spin band, and ii) it opens 
a conduction channel at the Fermi level in the minority 
spin band. This is a crucial aspect for understanding the 
transport properties of such structures. 
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FIG. 8: Ballistic conductance as a function of energy for a 50% 
Mn DFH with one AsGa in the middle of the unit cell, half the 
way between two consecutive (Ga,Mn)As planes. The vertical 
line indicates the position of the Fermi level, (a) majority, (b) 
minority spins. 



uniform system can be obtained simply by summing up 
the number of scattering channels at the Fermi level 

2 bz 2 BZ 2 

r 5 4E ] 4E N k< = t ( Nv ) a > (6) 



B. Transport in presence of As antisites 

We calculate the in plane ballistic conductance for a 
50% Mn DFH with one Asca in the middle of the spacer 
and the Mn ions uniformly distributed in the plane (con- 
figuration FMi). This is the same situation as in figure|S] 
The results are presented in figure^ where again we have 
considered up to 100 /c-points in the transverse Brillouin 
zone. 

The main difference with respect to the AsQa-free case 
is the presence of a non-vanishing conductance for the 
minority spin electrons at the Fermi energy. This is en- 
tirely due to the Asc a impurity band as pointed out in the 
previous section, and of course destroys the half-metallic 
behavior of the spin current. The spin-polarization of the 
current P, defined as 



P : 



(5) 



where the conductances are taken at the Fermi level, 
drops from 100% in the defect-free case to about 60% 
in the present case. 

However the situation is rather different if we consider 
diffusive transport 46]. This is due to the different way 
in which the density of states and the group velocity con- 
tribute to the conductance in the ballistic and diffusive 
limits. In fact, in the ballistic limit the conductance of a 



where and v% are respectively the density of state and 
the group velocity for a spin a scattering channel, and 
the sum is performed over the two dimensional Brillouin 
zone orthogonal to the transport direction. In contrast 
from the classical Boltzmann equation one finds that the 
diffusive spin conductance is proportional to 



(7) 



where now the sum runs over the three dimensional Fermi 
surface and t° is the spin-dependent relaxation time 0] . 
This reflects the well-known fact that, while in the ballis- 
tic limit all the scattering channels contribute with e 2 /h 
to the conductance independently from their group ve- 
locity, in the diffusive case the current is dominated by 
fast electrons because of the v 2 dependence. 

If one assumes that the relaxation time is not depen- 
dent on the spin direction (r^ = t^), then the spin polar- 
ization of the current in the diffusive limit can be written 



(Nv 2 ) 



(Nv 2 y - {Nv 2 )± 
(Nv 2 )1 + (Nv 2 )^ 



(8) 



In general P and P^u 2 ) are different, with P^u 2 ) larger 
if the difference in conductance for the two spin channels 
originates from a large Fermi velocity mismatch between 
the two spin bands. This is the case in the present DFH. 



In the majority band the Fermi surface is derived from 
the top of the GaAs valence band and the Fermi velocity 
is rather large. In contrast the minority band Fermi sur- 
face is due to the As antisite impurity band and the Fermi 
velocity is quite small. Therefore, since in the diffusive 
limit the conductance is proportional to v 2 , we expect 
a much larger spin polarization of the current compared 
with the ballistic case. 




FIG. 9: (Nv 2 ) as a function of energy for a 50% Mn DFH 
with one AsGa in the middle of the unit cell, half the way 
between two consecutive (Ga,Mn)As planes. The vertical line 
indicates the position of the Fermi level, (a) majority, (b) 
minority spins. 

In figure [§] we present (Nv 2 ) as a function of energy 
for the same DFH as that of figure |H1 The spin polar- 
ization at the Fermi energy is now 80%. Since in actual 
DFH the transport is due to hopping conductance [l6[ . 
we can conclude that As antisites, although they destroy 
the half-metallic state do not strongly affect the spin- 
polarization of the current. 

Finally it is interesting to study the spatial distribu- 
tion of the current across the DFH in the presence of As 
antisites. In figure ITU1 we present the real space charge 
density distribution, p(r), calculated only for those states 
contributing to the conductance at the Fermi energy • 
Strictly speaking this does not represent the current dis- 
tribution in real space, but gives information on the spa- 
tial distribution of the conductance electrons' wave func- 
tions around E-p . The main feature of figure 1101 is that 
there is a spacial separation between the two spin cur- 
rents, with the majority spin current located near the Mn 
plane, and the minority current strongly localized around 
the As antisites. The first is evenly distributed and this 
is suggestive of a metallic like behavior, while the second 
is strongly localized at the scattering center, suggesting 
an hopping-likc transport. Moreover, if we correlate the 
spin current with the relevant band-structure (see figure 
0), we notice that the majority spin-current is hole- like, 
while the minority is electron-like. 
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FIG. 10: Real space distribution of the current. This is calcu- 
lated as the charge density distribution in real space of those 
scattering states contributing to the conductance at Ef 



We also investigate how these features change when 
the As antisite moves toward the Mn plane. Although the 
conductance as a function of energy does not present any 
significant modifications with respect to the case of figure 
\7\ the spacial distribution shows an increased participa- 
tion of the As antisite in the majority spin conductance 
as it moves closer to the Mn plane. This suggests that 
As antisites can play an important role in spin relaxation 
processes within DFH. 

In conclusion, our transport results suggest that, if 
there is no spin mixing, the transport is dominated by a 
p-type metallic-like majority spin current with a smaller 
contribution from an n-type hopping-like minority spin 
current. Moreover the two spin-currents are spatially sep- 
arated, preventing spin-mixing, only if the As antisites 
are reasonably far from the Mn plane. 



C. Fit to the Heisenberg Model 

We now investigate the stability of the ferromagnetic 
state when As antisites are present. We perform similar 
calculations to those described in section IIIC, but this 
time we include one Asc a in the unit cell. Since the rel- 
ative position of the Asc a with respect to the Mn ions is 
crucial in determining the electronic and magnetic prop- 
erties |35| . we investigate how the coupling depends on 
the As antisite position along the superlattice direction. 

In figure ITTk we present the value of the exchange con- 
stant for first, second and third nearest neighbor cou- 
plings as a function of the position of the As antisite 
with respect to the Mn plane. In the same figure we also 
present the same values for the As antisite-free case and 
for the case of two As antisites (located respectively at 
1/3 and 2/3 of the supercell along the superlattice direc- 
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tion). 
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FIG. 11: a) Exchange constants Ji, J2 and J3 and b) mag- 
netic energy -E mag of a DFH as a function of the position of 
the As antisite with respect to the Mn plane (in unit of the 
superlattice period). We also include the case of two As an- 
tisites. The horizontal lines indicate the same quantities for 
the antisite-free case. 

From the figure it is clear that, in all cases, most of 
the coupling comes from nearest neighbor interaction, 
which accounts on average for 75% of the total coupling 
Jtot = J2i Ji- Then there is a fast decay of the exchange 
coupling with separation. This coupling is therefore short 
range. Furthermore it is very interesting to note that the 
second nearest neighbor coupling (Mn ions at Al and 
A2) is almost negligible in all cases, and it sometimes 
assumes negative values. This suggests possible local an- 
tiferromagnetic coupling between the Mn ions as recently 
demonstrated experimentally |l5j . 

Turning our attention to the effect of As antisites, it is 
clear that these weaken the coupling between the Mn 
ions. This is expected since an AsQa is a donor and 
therefore its presence partially compensates the free holes 
responsible for the long range ferromagnetic coupling. 
However our accurate supercell calculations bring addi- 
tional interesting features. In agreement with the case of 
the random alloys (see reference |35j), the actual position 
of the As antisite with respect to the Mn ions determines 
most of the properties. From figure fTTI it is clear that the 
ferromagnetic coupling is weakened when the As antisite 
moves toward the Mn plane. In particular the exchange 
constants for the cases when z = 1/6 and z = 0, are very 
similar to those of the case of two As antisites (z = 1/3 
and z = 2/3). This feature can be understood by con- 
sidering the potential profile discussed in the previous 
section. In fact the charge is strongly confined in the Mn 
plane with a potential barrier separating the plane from 
the GaAs region. Therefore it is natural to think that 
any perturbation in the GaAs region will only weakly af- 
fect the electronic configuration of the Mn plane, unless 



such a perturbation is spatially located in the vicinity of 
such a plane. 

In the second half of figure^] we present the magnetic 
energy, E mag — x^ ji J i , for the different cells investi- 
gated. As we pointed out previously this quantity is pro- 
portional to the Curie temperature, Tq- From the figure 
it is clear that E mag depends sensitively on the presence 
of the As antisites and on their actual location. In par- 
ticular, although we find the lowest value of E mag in the 
case in which two As antisites are present, we also find 
that this is very similar to the case of a single As antisite 
located in close proximity to the Mn plane [z = 1/6 in 
the present case). Therefore we conclude that one can 
obtain high Curie temperatures, not only by avoiding 
the formation of intrinsic defects, but also by controlling 
their position with respect to the magnetic region. 

Finally we make a few comments on the effect of disor- 
der. Since in actual samples we are not able to control the 
exact position of the Mn ions with respect to each other, 
and since the Mn-Mn coupling is strongly dependent on 
the relative positions of the Mn ions and those of the As 
antisites, it is likely that there are regions of strong Mn- 
Mn coupling together with regions of weak or even an- 
tiferromagnetic coupling. This suggests that in addition 
to configurational disorder, magnetic disorder can also 
be present in DFH even at low temperatures. Therefore 
since the electrons (or holes) at the Fermi level have a 
rather large (metallic) density, and are strongly confined 
in a few atomic planes around the Mn ions, we can con- 
clude that DFH have the electronic properties of highly 
resistive (dirty) metals. Turning the argument around, 
we conclude that the metallicity is crucial for the mag- 
netic state of DFH, and that the most metallic samples 
are likely to show less magnetic disorder and therefore 
more robust magnetic properties |17|. 



V. AlAs/(Ga,Mn)As HETEROSTRUCTURES 

One of the main messages from the analysis done so far 
is that the exchange part of the DFT potential creates a 
strong confinement potential for the majority electrons in 
the Mn plane. In this section we investigate the effects of 
an additional confining potential, namely that obtained 
by sandwiching a (Ga,Mn)As monolayer into the GaAs 
region of an AlAs/GaAs superlattice. Our expectation is 
that the AlAs/GaAs band alignment will further confine 
the spin-holes in the proximity of the Mn ions, therefore 
enhancing the exchange coupling. 



A. GaAs/AlAs band alignment 

Before considering the Mn-doped case we first illus- 
trate the general band alignment of an AlAs/GaAs su- 
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perlattice, as obtained from our DFT calculations. The 
valence band off-set, A, is calculated as suggested by Ba- 
roni et al. |47j as 



AE V + AV 



(9) 



where AV is the off-set between the AlAs and the GaAs 
electrostatic potentials calculated for an AlAs / GaAs het- 
erostructure, and AE V is the energy difference between 
the valence band tops, E v , of the bulk materials. These 
are calculated from their electrostatic potential V 



AE V = (E v - ^) GaAs - (E v - V) 



AlAs 



(10) 



Here we have constructed a supercell by stacking four 
AlAs cubic cells on top of four GaAs cubic cells, all with 
the same GaAs lattice constant (5.65A). Our calculations 
give a valence band off-set of A=0.405 eV (A=0.403 eV if 
the total DFT potential is considered instead of the elec- 
trostatic one) . This value for A is in very good agreement 
with both experimental results |48| and early DFT cal- 
culations [49(. The resulting GaAs/ AlAs band alignment 
is shown in figure 1121 it provides an additional confine- 
ment potential for holes in the GaAs region. Therefore a 
stronger ferromagnetic coupling is expected. 
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FIG. 12: GaAs/AlAs band alignment. 

B. Electronic Structure 

Also in this case the band structure and the DFT 
potential are the main quantities to investigate. First 
we study the evolution of the electronic properties as a 
function of the AlAs fraction of the superlattice. We con- 
struct GaAs(5_„)/ 2 /AlAs„/GaAs(5_„)/2/(Ga,Mn)Asi 
superlattices, where the labels indicate the number of 
monolayers of the specific semiconductor. Note that the 



total number of monolayers in our supercell is six, and 
that the AlAs fraction is always located in the middle of 
the cell. 
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FIG. 13: Band structure of AlAs/GaAs/(Ga,Mn)As DFH as 
a function of the AlAs fraction: a) n—0, b) n=l, c) n—3, d) 
n—5. The top panels are for the majority spin band and the 
bottom for the minority. The horizontal line indicates the 
position of the Fermi level, which is set to Ef=0 eV. 

In figure El we present the band structure of a 50% 
DFH with Mn ions uniformly distributed in the plane 
(configuration FMi), for the cases n=0,l,3,5. As usual 
we consider only the in plane directions. 

From figure ^] we can see that the only appreciable 
change in the band structure as the AlAs fraction in- 
creases is an enhancement of the band gap. This of course 
reflects the different AlAs/ GaAs ratios of the superlat- 
tices and the larger bandgap of AlAs compared with 
GaAs (our LDA values are 1.7 eV and 0.6 eV respec- 
tively). For all the superlattices studied the magnetic 
moment of the unit cell is 16 fj,B and the Miilliken analy- 
sis gives a Mn d occupation of ~4.7 and ~0.8 respectively 
for majority and minority spins. Such an occupation is 
independent of the AlAs fraction and suggests that Mn is 
always in a d 5 state plus an antiferromagnetically coupled 
hole, as in the case of GaAs/(Ga,Mn)As DFH. 

A closer look at the band structure reveals another 
important feature. The spin splitting of the valence band 
top, A„ — El — E~ , is a non-monotonic function of the 
AlAs fraction, with values of 0.88 eV, 0.93 eV, 1.02 eV, 
and 0.96 eV for n =0, 1, 3 and 5 respectively. The initial 
increase is due to the enhanced confinement of the free 
holes in the Mn region. In fact the mean field expression 
for the valence band top spin splitting is simply [J, |3Jj 



A„ = xN/3(S) 



(11) 



where N(3 is the p-d exchange constant, (S) is the mean 
spin ((S)=5/2 for (Ga,Mn)As) and x is the Mn concen- 
tration. The exchange constant, Nf3, depends on the 



12 



degree of overlap between the hole density and the Mn 
ions, and this is enhanced by hole confinement. There- 
fore we expect an increase of A„ when the AlAs fraction 
is increased. The case n = 5 is different, since no GaAs 
region is left and the additional confinement of the hole in 
the Mn plane due to the AlAs/GaAs valence band offset 
is partially lost. 

Another important quantity to investigate is the 
DFT total potential. In figure [21 we present the 
macroscopic average along the superlattice direction 
of the DFT and Hartree potentials, and the elec- 
tronic charge density, for a 50% DFH (configura- 
tion FMi), where we introduce three AlAs monolayers 
(GaAsi/AlAs 3 /GaAsi/(Ga,Mn)Asi). The figure shows 



1 

- . 1 




/^\'~y^\ — ' Spin 
\ — 4- Spin 

i,i, ,i 


i 1 i 
i , i 


(a): 


- 1 1 
1 >l 


1 > "N ' ' 1 
i , i \ 1 , / 


i , i 


(b): 


= • Mr 

- ♦ As 
lw Al 

i 1 * Ga 


1 r\ 1 7 V 

M M (C): 





4 8 12 16 20 24 28 32 



Position (A) 

FIG. 14: Macroscopic average of a) the total DFT potential, 
b) the Hartree potential and the c) charge density distribu- 
tion as a function of the position along the superlattice di- 
rection. The system is a GaAsi/AlAs3/GaAsi/(Ga,Mn)Asi 
superlattice with 50% Mn uniformely distributed in the plane 
(configuration FMi). Here we plot two superlattice periods. 



that there is a large well for both majority and minor- 
ity spin electrons in the AlAs region. This is mainly 
due to the Hartree component of the DFT potential and 
it is not spin sensitive. The width of this region de- 
pends on the thickness of the AlAs layer and plots for 
other values of n give similar results. If we interpret the 
macroscopic average of the DFT potential in the spirit 
of the envelope function approximation, we can conclude 
that new electrons coming from intrinsic donors will be 
preferentially localized in the AlAs region. Therefore 
the ferromagnetism in AlAs/(Ga,Mn)As DFH appears 
to be more robust against electron doping than that 
in GaAs/(Ga,Mn)As DFH. This is consistent with the 
rather large T c found in AlAs/GaAs/ (Ga,Mn) As DFH 

mm. 



C. Fit to the Heisenberg Model 

As in the case of (Ga,Mn) As / GaAs DFH we investigate 
the strength of the ferromagnetic coupling by fitting our 
total energy calculation to third nearest neighbor Heisen- 
berg model. We consider only the case of 50% Mn in the 
plane and we study the dependence of the coupling con- 
stants on the number of AlAs layers in the supercell. The 
results of our fit are shown in table Hill 

TABLE III: Ji , J2 , J3 and x ^% Ji as a function of the number 
of the AlAs layers in the supercell. All the units are meV. 



AlAs planes (n) 


Ji 


J 2 


■h 




n = 


13.3 


0.9 


4.5 


9.35 


n = 2 


13.8 


0.9 


4.7 


9.72 


n — 3 


14.0 


0.8 


4.9 


9.85 


n = 5 


14.6 


0.6 


5.2 


10.30 



From the table one can conclude that the total strength 
of the coupling, which is proportional to Tc, increases 
as a function of the AlAs fraction. This is consistent 
with an enhanced ferromagnetic coupling between the 
holes and the Mn ions due to confinement. Therefore in 
the absence of intrinsic defects, AlAs/(Ga,Mn)As DFH 
are more promising than GaAs/ (Ga,Mn)As DFH as high 
Tc DMS. It is also interesting to note that the different 
exchange constants behave in a different way when the 
AlAs fraction is increased: J\ and J3 are enhanced and 
J2 is reduced. This behavior is not surprising within a 
carrier mediated ferromagnetism model (e.g. RKKY), 
since an increase of the AlAs fraction changes the value 
of the Fermi wave vector, and therefore the period of the 
exchange coupling. 

Having established that defect- free AlAs/(Ga,Mn)As 
DFH present higher T c than their GaAs/(Ga,Mn)As 
counterparts, we finally investigate the stability of the 
ferromagnetic coupling against the presence of As anti- 
sites. As usual we introduce one Asc a per unit cell at 
various positions along the superlattice direction. Here 
we consider two limiting cases: 1) n — 3 and the As anti- 
site is located in one of the GaAs planes adjacent to the 
Mn plane, 2) n = 5 and the As antisite is located in the 
middle of the unit cell (in the AlAs region). 

In the first case, we find the coupling constants to 
be Ji = 7.0 meV, J 2 = 0.0 meV, J 3 = 2.9 meV, 
x Ji — 4.95 meV. This is a considerable reduction 
of the ferromagnetic coupling with respect to the defect- 
free case. In this case the values of the coupling constants 
are almost identical to those of GaAs/(Ga,Mn)As DFH 
in which an antisite is introduced at the same position 
(in figure ITU the Asca position is 1/6). Therefore, if the 
As antisites are introduced in the GaAs region, there will 
be no beneficial effects from the AlAs/GaAs band align- 
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ment. 

The situation is rather different in the second case 
where the As antisites are introduced in the AlAs layer. 
Now the coupling constants are J\ — 10.1 meV, J2 = 
—0.1 meV, J3 = 4.4 meV, and the reduction of the to- 
tal coupling x J2i Ji with respect to the defect-free case 
is only of about 30%. Interestingly J 2 now assumes a 
negative value, suggesting some possible frustration even 
at low temperature. 

In conclusion, AlAs/GaAs/(Ga,Mn)As DFH have 
stronger ferromagnetic interaction between the Mn ions 
than GaAs/(Ga,Mn)As DFH. Moreover, in the case that 
intrinsic defects are kept into the AlAs region, the ferro- 
magnetic order is more robust against hole compensation. 

VI. CONCLUSIONS 

We have performed an extensive theoretical study 
of the electronic, magnetic and transport properties 
of GaAs / (Ga,Mn) As and AlAs/GaAs/(Ga,Mn)As DFH, 
using DFT within LSDA. 

We find that GaAs/ (Ga,Mn)As DFH show an half- 
metal band structure with metallic conductance in the 
Mn plane. The macroscopic average of the DFT potential 
indicates a selective confinement of the spin holes in the 
Mn planes and the Mn-Mn ferromagnetic interaction is 
non-monotonically dependent on the Mn concentration. 

When compensating defects such as As antisites are 
introduced, the half-metallic state is lost and conduct- 
ing channels appear in the minority spin band. These 
are due to hopping conductance through localized AsQa 
states. However, at least when the As antisites are far 
from the Mn planes, there is a spatial separation of the 
two spin currents with a metallic majority spin current 



located in the Mn planes, and an hopping-type minority 
spin current located primarily in the GaAs region. These 
differences in the type of transport for the two spin bands 
are magnified in the diffusive limit, for which we calculate 
a spin polarization of about 80%. Finally, the presence 
of As antisites generally weakens the ferromagnetic in- 
teraction, and local antiferromagnetic coupling between 
Mn ions is possible at low temperature. 

With all these results in hand we conclude that 
GaAs/(Ga,Mn)As DFH behave as dirty planar metals, 
where the strength of the ferromagnetic coupling depends 
strongly on the amount and the position of the intrinsic 
defects. 

Finally we have investigated the effect of additional 
confinement by studying AlAs/GaAs/(Ga,Mn)As DFH. 
In this case our accurate total energy calculations confirm 
that the band offset between GaAs and AlAs strongly 
confines the holes in the Mn region resulting in a larger 
Mn-Mn coupling. In addition we find that these struc- 
tures are less sensitive to the presence of As antisites than 
the simpler GaAs/ (Ga,Mn)As DFH. 

We therefore conclude that band engineering, selective 
doping and the ability to control the actual position of As 
antisites or other compensating defects are the main keys 
to obtain room temperature ferromagnetism in DFH. 
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